Marine gregarine genomes reveal the breadth of apicomplexan diversity with a partially conserved glideosome machinery.

Our current view of the evolutionary history, coding and adaptive capacities of Apicomplexa, protozoan parasites of a wide range of metazoan, is currently strongly biased toward species infecting humans, as data on early diverging apicomplexan lineages infecting invertebrates is extremely limited. Here, we characterized the genome of the marine eugregarine Porospora gigantea, intestinal parasite of Lobsters, remarkable for the macroscopic size of its vegetative feeding forms (trophozoites) and its gliding speed, the fastest so far recorded for Apicomplexa. Two highly syntenic genomes named A and B were assembled. Similar in size (~ 9 Mb) and coding capacity (~ 5300 genes), A and B genomes are 10.8% divergent at the nucleotide level, corresponding to 16-38 My in divergent time. Orthogroup analysis across 25 (proto)Apicomplexa species, including Gregarina niphandrodes, showed that A and B are highly divergent from all other known apicomplexan species, revealing an unexpected breadth of diversity. Phylogenetically these two species branch sisters to Cephaloidophoroidea, and thus expand the known crustacean gregarine superfamily. The genomes were mined for genes encoding proteins necessary for gliding, a key feature of apicomplexans parasites, currently studied through the molecular model called glideosome. Sequence analysis shows that actin-related proteins and regulatory factors are strongly conserved within apicomplexans. In contrast, the predicted protein sequences of core glideosome proteins and adhesion proteins are highly variable among apicomplexan lineages, especially in gregarines. These results confirm the importance of studying gregarines to widen our biological and evolutionary view of apicomplexan species diversity, and to deepen our understanding of the molecular bases of key functions such as gliding, well known to allow access to the intracellular parasitic lifestyle in Apicomplexa.

Intracellular development and impact of a marine eukaryotic parasite on its zombified microalgal host.

Parasites are widespread and diverse in oceanic plankton and many of them infect single-celled algae for survival. How these parasites develop and scavenge energy within the host and how the cellular organization and metabolism of the host is altered remain open questions. Combining quantitative structural and chemical imaging with time-resolved transcriptomics, we unveil dramatic morphological and metabolic changes of the marine parasite Amoebophrya (Syndiniales) during intracellular infection, particularly following engulfment and digestion of nutrient-rich host chromosomes. Changes include a sequential acristate and cristate mitochondrion with a 200-fold increase in volume, a 13-fold increase in nucleus volume, development of Golgi apparatus and a metabolic switch from glycolysis (within the host) to TCA (free-living dinospore). Similar changes are seen in apicomplexan parasites, thus underlining convergent traits driven by metabolic constraints and the infection cycle. In the algal host, energy-producing organelles (plastid, mitochondria) remain relatively intact during most of the infection. We also observed that sugar reserves diminish while lipid droplets increase. Rapid infection of the host nucleus could be a "zombifying" strategy, allowing the parasite to digest nutrient-rich chromosomes and escape cytoplasmic defense, whilst benefiting from maintained carbon-energy production of the host cell.

Parasitic infections by Group II Syndiniales target selected dinoflagellate host populations within diverse protist assemblages in a model coastal pond.

Protists are integral to marine food webs and biogeochemical cycles; however, there is a paucity of data describing specific ecological niches for some of the most abundant taxa in marker gene libraries. Syndiniales are one such group, often representing the majority of sequence reads recovered from picoplankton samples across the global ocean. However, the prevalence and impacts of syndinian parasitism in marine environments remain unclear. We began to address these critical knowledge gaps by generating a high-resolution time series (March-October 2018) in a productive coastal pond. Seasonal shifts in protist populations, including parasitic Syndiniales, were documented during periods of higher primary productivity and increased summer temperature-driven stratification. Elevated concentrations of infected hosts and free-living parasite spores occurred at nearly monthly intervals in July, August, and September. We suggest intensifying stratification during this period correlated with the increased prevalence of dinoflagellates that were parasitized by Group II Syndiniales. Infections in some protist populations were comparable to previously reported large single-taxon dinoflagellate blooms. Infection dynamics in Salt Pond demonstrated the propagation of syndinian parasites through mixed protist assemblages and highlighted patterns of host/parasite interactions that better reflect many other marine environments where single taxon blooms are uncommon.

Temperature Affects the Biological Control of Dinoflagellates by the Generalist Parasitoid Parvilucifera rostrata.

The increase in emerging harmful algal blooms in the last decades has led to an extensive concern in understanding the mechanisms behind these events. In this paper, we assessed the growth of two blooming dinoflagellates (Alexandrium minutum and Heterocapsa triquetra) and their susceptibility to infection by the generalist parasitoid Parvilucifera rostrata under a temperature gradient. The growth of the two dinoflagellates differed across a range of temperatures representative of the Penzé Estuary (13 to 22 °C) in early summer. A. minutum growth increased across this range and was the highest at 19 and 22 °C, whereas H. triquetra growth was maximal at intermediate temperatures (15-18 °C). Interestingly, the effect of temperature on the parasitoid infectivity changed depending on which host dinoflagellate was infected with the dinoflagellate responses to temperature following a positive trend in A. minutum (higher infections at 20-22 °C) and a unimodal trend in H. triquetra (higher infections at 18 °C). Low temperatures negatively affected parasitoid infections in both hosts (i.e., "thermal refuge"). These results demonstrate how temperature shifts may not only affect bloom development in microalgal species but also their control by parasitoids.

From the sxtA4 Gene to Saxitoxin Production: What Controls the Variability Among Alexandrium minutum and Alexandrium pacificum Strains?

Paralytic shellfish poisoning (PSP) is a human foodborne syndrome caused by the consumption of shellfish that accumulate paralytic shellfish toxins (PSTs, saxitoxin group). In PST-producing dinoflagellates such as Alexandrium spp., toxin synthesis is encoded in the nuclear genome via a gene cluster (sxt). Toxin production is supposedly associated with the presence of a 4th domain in the sxtA gene (sxtA4), one of the core genes of the PST gene cluster. It is postulated that gene expression in dinoflagellates is partially constitutive, with both transcriptional and post-transcriptional processes potentially co-occurring. Therefore, gene structure and expression mode are two important features to explore in order to fully understand toxin production processes in dinoflagellates. In this study, we determined the intracellular toxin contents of twenty European Alexandrium minutum and Alexandrium pacificum strains that we compared with their genome size and sxtA4 gene copy numbers. We observed a significant correlation between the sxtA4 gene copy number and toxin content, as well as a moderate positive correlation between the sxtA4 gene copy number and genome size. The 18 toxic strains had several sxtA4 gene copies (9-187), whereas only one copy was found in the two observed non-toxin producing strains. Exploration of allelic frequencies and expression of sxtA4 mRNA in 11 A. minutum strains showed both a differential expression and specific allelic forms in the non-toxic strains compared with the toxic ones. Also, the toxic strains exhibited a polymorphic sxtA4 mRNA sequence between strains and between gene copies within strains. Finally, our study supported the hypothesis of a genetic determinism of toxin synthesis (i.e., the existence of several genetic isoforms of the sxtA4 gene and their copy numbers), and was also consistent with the hypothesis that constitutive gene expression and moderation by transcriptional and post-transcriptional regulation mechanisms are the cause of the observed variability in the production of toxins by A. minutum.

Rapid protein evolution, organellar reductions, and invasive intronic elements in the marine aerobic parasite dinoflagellate Amoebophrya spp.

BACKGROUND: Dinoflagellates are aquatic protists particularly widespread in the oceans worldwide. Some are responsible for toxic blooms while others live in symbiotic relationships, either as mutualistic symbionts in corals or as parasites infecting other protists and animals. Dinoflagellates harbor atypically large genomes (~ 3 to 250 Gb), with gene organization and gene expression patterns very different from closely related apicomplexan parasites. Here we sequenced and analyzed the genomes of two early-diverging and co-occurring parasitic dinoflagellate Amoebophrya strains, to shed light on the emergence of such atypical genomic features, dinoflagellate evolution, and host specialization. RESULTS: We sequenced, assembled, and annotated high-quality genomes for two Amoebophrya strains (A25 and A120), using a combination of Illumina paired-end short-read and Oxford Nanopore Technology (ONT) MinION long-read sequencing approaches. We found a small number of transposable elements, along with short introns and intergenic regions, and a limited number of gene families, together contribute to the compactness of the Amoebophrya genomes, a feature potentially linked with parasitism. While the majority of Amoebophrya proteins (63.7% of A25 and 59.3% of A120) had no functional assignment, we found many orthologs shared with Dinophyceae. Our analyses revealed a strong tendency for genes encoded by unidirectional clusters and high levels of synteny conservation between the two genomes despite low interspecific protein sequence similarity, suggesting rapid protein evolution. Most strikingly, we identified a large portion of non-canonical introns, including repeated introns, displaying a broad variability of associated splicing motifs never observed among eukaryotes. Those introner elements appear to have the capacity to spread over their respective genomes in a manner similar to transposable elements. Finally, we confirmed the reduction of organelles observed in Amoebophrya spp., i.e., loss of the plastid, potential loss of a mitochondrial genome and functions. CONCLUSION: These results expand the range of atypical genome features found in basal dinoflagellates and raise questions regarding speciation and the evolutionary mechanisms at play while parastitism was selected for in this particular unicellular lineage.

Dinophyceae can use exudates as weapons against the parasite Amoebophrya sp. (Syndiniales).

Parasites in the genus Amoebophrya sp. infest dinoflagellate hosts in marine ecosystems and can be determining factors in the demise of blooms, including toxic red tides. These parasitic protists, however, rarely cause the total collapse of Dinophyceae blooms. Experimental addition of parasite-resistant Dinophyceae (Alexandrium minutum or Scrippsiella donghaienis) or exudates into a well-established host-parasite coculture (Scrippsiella acuminata-Amoebophrya sp.) mitigated parasite success and increased the survival of the sensitive host. This effect was mediated by waterborne molecules without the need for a physical contact. The strength of the parasite defenses varied between dinoflagellate species, and strains of A. minutum and was enhanced with increasing resistant host cell concentrations. The addition of resistant strains or exudates never prevented the parasite transmission entirely. Survival time of Amoebophrya sp. free-living stages (dinospores) decreased in presence of A. minutum but not of S. donghaienis. Parasite progeny drastically decreased with both species. Integrity of the dinospore membrane was altered by A. minutum, providing a first indication on the mode of action of anti-parasitic molecules. These results demonstrate that extracellular defenses can be an effective strategy against parasites that protects not only the resistant cells producing them, but also the surrounding community.

Unveiling protist diversity associated with the Pacific oyster Crassostrea gigas using blocking and excluding primers.

BACKGROUND: Microbiome of macroorganisms might directly or indirectly influence host development and homeostasis. Many studies focused on the diversity and distribution of prokaryotes within these assemblages, but the eukaryotic microbial compartment remains underexplored so far. RESULTS: To tackle this issue, we compared blocking and excluding primers to analyze microeukaryotic communities associated with Crassostrea gigas oysters. High-throughput sequencing of 18S rRNA genes variable loops revealed that excluding primers performed better by not amplifying oyster DNA, whereas the blocking primer did not totally prevent host contaminations. However, blocking and excluding primers showed similar pattern of alpha and beta diversities when protist communities were sequenced using metabarcoding. Alveolata, Stramenopiles and Archaeplastida were the main protist phyla associated with oysters. In particular, Codonellopsis, Cyclotella, Gymnodinium, Polarella, Trichodina, and Woloszynskia were the dominant genera. The potential pathogen Alexandrium was also found in high abundances within some samples. CONCLUSIONS: Our study revealed the main protist taxa within oysters as well as the occurrence of potential oyster pathogens. These new primer sets are promising tools to better understand oyster homeostasis and disease development, such as the Pacific Oyster Mortality Syndrome (POMS) targeting juveniles.

Identification to species level of live single microalgal cells from plankton samples with matrix-free laser/desorption ionization mass spectrometry.

INTRODUCTION: Marine planktonic communities are complex microbial consortia often dominated by microscopic algae. The taxonomic identification of individual phytoplankton cells usually relies on their morphology and demands expert knowledge. Recently, a live single-cell mass spectrometry (LSC-MS) pipeline was developed to generate metabolic profiles of microalgae. OBJECTIVE: Taxonomic identification of diverse microalgal single cells from collection strains and plankton samples based on the metabolic fingerprints analyzed with matrix-free laser desorption/ionization high-resolution mass spectrometry. METHODS: Matrix-free atmospheric pressure laser-desorption ionization mass spectrometry was performed to acquire single-cell mass spectra from collection strains and prior identified environmental isolates. The computational identification of microalgal species was performed by spectral pattern matching (SPM). Three similarity scores and a bootstrap-derived confidence score were evaluated in terms of their classification performance. The effects of high and low-mass resolutions on the classification success were evaluated. RESULTS: Several hundred single-cell mass spectra from nine genera and nine species of marine microalgae were obtained. SPM enabled the identification of single cells at the genus and species level with high accuracies. The receiver operating characteristic (ROC) curves indicated a good performance of the similarity measures but were outperformed by the bootstrap-derived confidence scores. CONCLUSION: This is the first study to solve taxonomic identification of microalgae based on the metabolic fingerprints of the individual cell using an SPM approach.

Cryptic species in the parasitic Amoebophrya species complex revealed by a polyphasic approach.

As critical primary producers and recyclers of organic matter, the diversity of marine protists has been extensively explored by high-throughput barcode sequencing. However, classification of short metabarcoding sequences into traditional taxonomic units is not trivial, especially for lineages mainly known by their genetic fingerprints. This is the case for the widespread Amoebophrya ceratii species complex, parasites of their dinoflagellate congeners. We used genetic and phenotypic characters, applied to 119 Amoebophrya individuals sampled from the same geographic area, to construct practical guidelines for species delineation that could be applied in DNA/RNA based diversity analyses. Based on the internal transcribed spacer (ITS) regions, ITS2 compensatory base changes (CBC) and genome k-mer comparisons, we unambiguously defined eight cryptic species among closely related ribotypes that differed by less than 97% sequence identity in their SSU rDNA. We then followed the genetic signatures of these parasitic species during a three-year survey of Alexandrium minutum blooms. We showed that these cryptic Amoebophrya species co-occurred and shared the same ecological niche. We also observed a maximal ecological fitness for parasites having narrow to intermediate host ranges, reflecting a high cost for infecting a broader host range. This study suggests that a complete taxonomic revision of these parasitic dinoflagellates is long overdue to understand their diversity and ecological role in the marine plankton.

Dinoflagellate Host Chloroplasts and Mitochondria Remain Functional During Amoebophrya Infection.

Dinoflagellates are major components of phytoplankton that play critical roles in many microbial food webs, many of them being hosts of countless intracellular parasites. The phototrophic dinoflagellate Scrippsiella acuminata (Dinophyceae) can be infected by the microeukaryotic parasitoids Amoebophrya spp. (Syndiniales), some of which primarily target and digest the host nucleus. Early digestion of the nucleus at the beginning of the infection is expected to greatly impact the host metabolism, inducing the knockout of the organellar machineries that highly depend upon nuclear gene expression, such as the mitochondrial OXPHOS pathway and the plastid photosynthetic carbon fixation. However, previous studies have reported that chloroplasts remain functional in swimming host cells infected by Amoebophrya. We report here a multi-approach monitoring study of S. acuminata organelles over a complete infection cycle by nucleus-targeting Amoebophrya sp. strain A120. Our results show sustained and efficient photosystem II activity as a hallmark of functional chloroplast throughout the infection period despite the complete digestion of the host nucleus. We also report the importance played by light on parasite production, i.e., the amount of host biomass converted to parasite infective propagules. Using a differential gene expression analysis, we observed an apparent increase of all 3 mitochondrial and 9 out of the 11 plastidial genes involved in the electron transport chains (ETC) of the respiration pathways during the first stages of the infection. The longer resilience of organellar genes compared to those encoded by the nucleus suggests that both mitochondria and chloroplasts remain functional throughout most of the infection. This extended organelle functionality, along with higher parasite production under light conditions, suggests that host bioenergetic organelles likely benefit the parasite Amoebophrya sp. A120 and improve its fitness during the intracellular infective stage.

An aerobic eukaryotic parasite with functional mitochondria that likely lacks a mitochondrial genome.

Dinoflagellates are microbial eukaryotes that have exceptionally large nuclear genomes; however, their organelle genomes are small and fragmented and contain fewer genes than those of other eukaryotes. The genus Amoebophrya (Syndiniales) comprises endoparasites with high genetic diversity that can infect other dinoflagellates, such as those forming harmful algal blooms (e.g., Alexandrium). We sequenced the genome (~100 Mb) of Amoebophrya ceratii to investigate the early evolution of genomic characters in dinoflagellates. The A. ceratii genome encodes almost all essential biosynthetic pathways for self-sustaining cellular metabolism, suggesting a limited dependency on its host. Although dinoflagellates are thought to have descended from a photosynthetic ancestor, A. ceratii appears to have completely lost its plastid and nearly all genes of plastid origin. Functional mitochondria persist in all life stages of A. ceratii, but we found no evidence for the presence of a mitochondrial genome. Instead, all mitochondrial proteins appear to be lost or encoded in the A. ceratii nucleus.

First report of vampyrellid predator-prey dynamics in a marine system.

We report for the first time the in situ dynamics of a vampyrellid in a marine system. A high sampling frequency (twice-weekly) was applied in a tropical eutrophic lagoon (Rio de Janeiro, Brazil) for 5 years (2012-2016). The vampyrellid Hyalodiscus sp. specifically fed on the diatom Chaetoceros minimus during a short time window (~3 months), although the prey was intermittently detected as the dominant phytoplanktonic species over a longer period (~1 year). A classic Lotka-Volterra predator-prey dynamic was observed between the two partners, with a significant modification of the short-term oscillations of the prey. Specific abiotic preferences (i.e., relatively low temperature, intermediate salinity, and stratified conditions) associated with prey availability seemed to define this narrow temporal window of occurrence. Our results suggest that vampyrellids can be ecologically relevant in marine pelagic systems, with their impact on planktonic dynamics strongly depending on complex interactions between both biotic and abiotic factors.

Revisions to the Classification, Nomenclature, and Diversity of Eukaryotes.

This revision of the classification of eukaryotes follows that of Adl et al., 2012 [J. Euk. Microbiol. 59(5)] and retains an emphasis on protists. Changes since have improved the resolution of many nodes in phylogenetic analyses. For some clades even families are being clearly resolved. As we had predicted, environmental sampling in the intervening years has massively increased the genetic information at hand. Consequently, we have discovered novel clades, exciting new genera and uncovered a massive species level diversity beyond the morphological species descriptions. Several clades known from environmental samples only have now found their home. Sampling soils, deeper marine waters and the deep sea will continue to fill us with surprises. The main changes in this revision are the confirmation that eukaryotes form at least two domains, the loss of monophyly in the Excavata, robust support for the Haptista and Cryptista. We provide suggested primer sets for DNA sequences from environmental samples that are effective for each clade. We have provided a guide to trophic functional guilds in an appendix, to facilitate the interpretation of environmental samples, and a standardized taxonomic guide for East Asian users.

Fine structure and Molecular Phylogenetic Position of Two Marine Gregarines, Selenidium pygospionis sp. n. and S. pherusae sp. n., with Notes on the Phylogeny of Archigregarinida (Apicomplexa).

Archigregarines are a key group for understanding the early evolution of Apicomplexa. Here we report morphological, ultrastructural, and molecular phylogenetic evidence from two archigregarine species: Selenidium pygospionis sp. n. and S. pherusae sp. n. They exhibited typical features of archigregarines. Additionally, an axial row of vacuoles of a presumably nutrient distribution system was revealed in S. pygospionis. Intracellular stages of S. pygospionis found in the host intestinal epithelium may point to the initial intracellular localization in the course of parasite development. Available archigregarine SSU (18S) rDNA sequences formed four major lineages fitting the taxonomical affiliations of their hosts, but not the morphological or biological features used for the taxonomical revision by Levine (1971). Consequently, the genus Selenidioides Levine, 1971 should be abolished. The branching order of these lineages was unresolved; topology tests rejected neither para- nor monophyly of archigregarines. We provided phylogenies based on LSU (28S) rDNA and near-complete ribosomal operon (concatenated SSU, 5.8S, LSU rDNAs) sequences including S. pygospionis sequences. Although being preliminary, they nevertheless revealed the monophyly of gregarines previously challenged by many molecular phylogenetic studies. Despite their molecular-phylogenetic heterogeneity, archigregarines exhibit an extremely conservative plesiomorphic structure; their ultrastructural key features appear to be symplesiomorphies rather than synapomorphies.

Comparative Time-Scale Gene Expression Analysis Highlights the Infection Processes of Two Amoebophrya Strains.

Understanding factors that generate, maintain, and constrain host-parasite associations is of major interest to biologists. Although little studied, many extremely virulent micro-eukaryotic parasites infecting microalgae have been reported in the marine plankton. This is the case for Amoebophrya, a diverse and highly widespread group of Syndiniales infecting and potentially controlling dinoflagellate populations. Here, we analyzed the time-scale gene expression of a complete infection cycle of two Amoebophrya strains infecting the same host (the dinoflagellate Scrippsiella acuminata), but diverging by their host range (one infecting a single host, the other infecting more than one species). Over two-thirds of genes showed two-fold differences in expression between at least two sampled stages of the Amoebophrya life cycle. Genes related to carbohydrate metabolism as well as signaling pathways involving proteases and transporters were overexpressed during the free-living stage of the parasitoid. Once inside the host, all genes related to transcription and translation pathways were actively expressed, suggesting the rapid and extensive protein translation needed following host-cell invasion. Finally, genes related to cellular division and components of the flagellum organization were overexpressed during the sporont stage. In order to gain a deeper understanding of the biological basis of the host-parasitoid interaction, we screened proteins involved in host-cell recognition, invasion, and protection against host-defense identified in model apicomplexan parasites. Very few of the genes encoding critical components of the parasitic lifestyle of apicomplexans could be unambiguously identified as highly expressed in Amoebophrya. Genes related to the oxidative stress response were identified as highly expressed in both parasitoid strains. Among them, the correlated expression of superoxide dismutase/ascorbate peroxidase in the specialist parasite was consistent with previous studies on Perkinsus marinus defense. However, this defense process could not be identified in the generalist Amoebophrya strain, suggesting the establishment of different strategies for parasite protection related to host specificity.